function F = utility(x,gamma);

if gamma == 1
     F=log(x);
else
     F=((x)^(1-gamma)-1)/(1-gamma);    
end
end

